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We study experimentally and theoretically, the reorganization of superconducting vortices driven 
by oscillatory forces near the plastic depinning transition. We show that the system can be taken to 
configurations that are tagged by the shaking parameters but keep no trace of the initial conditions. 
In experiments performed in NbSe2 crystals, the periodic drive is induced by ac magnetic shaking 
fields and the overall order of the resulting configuration is determined by non invasive ac susceptibil- 
ity measurements. With a model of interacting particles driven over random landscapes, we perform 
molecular dynamics simulations that reveal the nature of the shaking dynamics as fluctuating states 
similar to those predicted for other interacting particle systems. 



PACS numbers: 74.25.Uv, 74.25.Wx, 74.62.En 

Vortices in type II superconductors 1 ^.!!, charge den- 
sity wavesi£ colloidal particles! 3 - - — , Wigner crystals** 
and domain walls!£ are only few of the many examples 
of physical systems belonging to the category of driven 
elastic manifolds over random landscapes. In all these 
systems, competing interactions give rise to a complex 
dynamics characterized by a region of plastic motion: un- 
der the action of an external driving force, strong enough 
to depin the system, a non-linear response develops, some 
particles are mobile while others remain pinned. This 
plastic depinning and its relation with the proliferation 
of topological defects has been object of study for several 
years. Recently, it has been claimed!^ that the plastic 
depinning of a colloidal system with random quenched 
disorder under the action of a dc force can be described 
in terms of a non equilibrium phase transition with a 
divergent transient time and it can be related to the 
"absorbing transition" from a random self-organization 
to a drive-dependent fluctuating steady state (FSS) ob- 
served in experiments on periodically sheared particle 
suspensions!^. 

Vortex matter offers an ideal playground to study 
this phenomenology both numerically and experimen- 
tally. Dynamic phase transition 1 ' 4 and plastic flow have 
been theoretically described and measured in the neigh- 
borhood of the order-disorder (O-D) transition, where 
the system evolves from a quasi-ordered Bragg glass 
(BG) 3 to a disordered phase with proliferation of topo- 
logical defects. The experimental fingerprint of the O-D 
transition is the anomalous non-monotonous dependence 
of the critical current density J c with both temperature 
and magnetic field, known as Peak Effect (PE)-. There, 
the response at a fixed temperature and field can de- 
pend on thermal, magnetic and dynamical history that 
can modify the topology of the vortex lattice configura- 
tions (VLCs) by creation or annihilation of disclinationsS.. 
In transport experiments vortex instability phenomena 
and the smearing of the PE have been claimed to orig- 
inate on a competition between the injection of a disor- 
dered vortex phase at the sample edges, and the dynamic 



annealing of this metastable disorder by the transport 
current^. However, using an ultrafast transport tech- 
nique that avoids current-induced vortex lattice (VL) re- 
organization, Xiao et al& have shown the existence of 
an enlarged crossing phase boundary between ordered 
and disordered phases in NbSe2 crystals arising from the 
bulk VL response, that has been corroborated in a re- 
cent work by our group using non invasive ac susceptibil- 
ity measurements in the linear regime^. In experiments, 
a shaking ac field is often applied to order the VL and 
the most accepted picture is that the shaking field assists 
the system in an equilibration process, from a disordered 
metastable configuration to the equilibrium BG phase. 
However, in our experiments, we have shown that in the 
transition region a shaking ac field can either order or 
disorder the VL, driving the system to VLCs with inter- 
mediate degree of disorder that are independent of the 
initial VLC. Very recently, noise transport experiments 
(in a-Mo^Gei-a;) in this intermediate region have shown 
evidence of a depinning transition with critical behavior 
similar to the absorbing transition 2 . 

In this work we study experimentally and numerically 
the reorganization of a vortex system in the plastic flow 
region under the action of ac shaking magnetic fields. We 
have done experiments changing the frequency and the 
waveform of the shaking fields and measuring the ac re- 
sponse of the final VLC in each case with very low ac 
amplitude, avoiding modifications in the final VLC. We 
have also performed numerical calculations that mimic 
our experiment. We have first used the Brandt method 18 
to determine the induced current dependence on the ex- 
ternal magnetic field and as a second step we have used 
this result as input for a molecular dynamics simulation. 
Although the model is rather simple, we obtain results 
which are qualitatively consistent with the experimen- 
tal results and show evidence of ac driven vortex lattice 
reorganization. 

Experimental results shown in this work correspond 
to a NbSe2 single crysta l 19 i 20 of approximate dimensions 
(0.5 x 0.5 x 0.03) mm 3 , with Tc = 7.30 K (defined as 
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FIG. 1. (Inset) Typical x" and x FCC and FCW curves in 
the linear regime measured at / = 30 kHz. Arrows indicate 
the direction of temperature variation. (Main panel) Zoom 
of the region delimited in the inset. Gray triangles (black 
circles) correspond to shaking frequency f s h = 100 kHz 
(fsh = 10 Hz). Open (full) symbols for warming (cooling) 
experiment. In darkgray small circles a FCC curve is shown. 
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FIG. 2. (a) x' versus time before and after shaking the sys- 
tem (the shaking time is not in scale and it is indicated by 
a gap without data). Gray (black) symbol correspond to 
f s h ~ 100 kHz (f s h = 10 Hz ). (b) x' values versus the 
shaking frequency (log scale). Arrows point from the values 
of x' before shaking to the values obtained after shaking, (c) 
X' values versus the fundamental shaking frequency for sinu- 
soidal (open dots) and square (black squares) shaking wave- 
forms. 



the mid point of the ac susceptibility linear transition 
at H = 0) and ATc = ±0.02 if. A neutron study^i 
in samples synthesized with the same technique as ours, 
shows an excellent agreement with the VL structure pre- 
dicted for randomly distributed point defects^. The ac 
susceptibility has been measured with a home made sus- 
ceptometer based in the mutual inductance technique, 
where both the ac and dc fields are parallel to the c axis 
of the sample. To study the VL response without dis- 
turbances, the measurements have been restricted to the 
linear Campbell regime^, where a very small ac field h a 
superimposed to the dc field H is applied, forcing vortices 
to perform small (harmonic) oscillations inside their ef- 
fective pinning potential wells. In this regime, the ac 
susceptibility is independent of amplitude and frequency 
and the inductive component of the ac susceptibility x' is 
determined by the experimental geometry and the curva- 
ture of the effective pinning potential well a l and a lower 
susceptibility x' (closer to -1) can be directly associated 
to a more pinned VL (see Ref. l24i ). Therefore, at fixed 
T and H, different x' (i- e - different degree of effective 
pinning) can be associated to the existence of VLCs with 
different degree of disorder. 

Figure Q] shows x' (T) measured in the linear regime 
(h a = 0.025 Oe,/ = 30 kHz), in a dc field H = 320 Oe 
after different thermal and dynamical histories. In the 
Inset, both components x' an d X" m field cool cool- 
ing (FCC) and field cool warming (FCW) processes are 
shown in the full temperature range, and the region of 
the PE (main panel) is indicated. The curve with small 
darkgray circles in the main panel was obtained in a FCC 
process, where the VL remains trapped in a metastable 



disordered and strong pinned configuration, and the PE 
disappears. The other experimental points have been 
obtained by measuring the system after shaking the VL 
at different frequencies. In all those cases, the follow- 
ing procedure has been performed: The system has been 
stabilized at a shaking temperature T s h and measured 
in the linear regime, recording the x' value correspond- 
ing to the starting configuration. Then the measurement 
was interrupted and N s h cycles of a sinusoidal shaking 
ac field (h s h = 3.2 Oe) were applied. After switching off 
the shaking field, x' corresponding to the final value was 
recorded in the linear regime; then the system was driven 
to the next shaking temperature. Full and open sym- 
bols correspond to procedures increasing and decreasing 
T s h- For clarity, in this figure only the final values are 
shown. Gray triangles indicate the response after shak- 
ing the system with N s h = 100 cycles at a frequency 
f s h = 100 kHz, whereas black circles show the response 
after a similar procedure at f s h, = 10 Hz. One remark- 
able feature published in our previous work4 is that the 
final VLCs do not show appreciable relaxation (within 
our experimental time window) and are independent on 
the initial condition. Surprisingly, we find that these fi- 
nal states are not unique: they depend on the shaking 
frequency, but however the independence of initial con- 
ditions persists. 

In Fig. d{a), the response at fixed temperature, before 
and after applying shaking protocols, is shown for differ- 
ent starting VLCs (same frequencies as in Fig. [I}. It can 
be seen that the shaking field can order or disorder the 
VL. The various curves X'(t) obtained at different initial 
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configurations collapse in an unique curve after shaking 
the system at a given frequency. Figure[3{b) shows the re- 
sponse after applying sinusoidal shaking fields at various 
frequencies. Arrows point from each starting VLC to the 
corresponding final configuration after shaking. In all the 
cases the final response is independent of the initial one, 
for shaking amplitudes greater than 1.8 Oe and above a 
number of applied cycles (in our system N s h ~ 40). For 
all the tested frequencies (not all shown in the figure), 
the higher the shaking frequency, the more ordered the 
resulting VLC, consistent with a dynamic reordering of 
the VLii. 

In that scenario, the frequency dependence would be 
due mainly to the different vortex velocities (i. e. forces) 
induced by the different shaking ac fields. Additional 
evidence supporting this fact arises from the response 
of the system to shaking fields with a square waveform. 
In this case it is expected that vortices move at higher 
velocities during short time intervals triggered and de- 
termined by the rising and falling edges, whereas the 
shaking frequency only modifies the waiting time be- 
tween these short intervals. In Fig^c), it is seen that for 
square shaking fields, the VLCs are maximally ordered. 
On the other hand, frequency independence shows that 
relaxation processes are not relevant at least in this time 
window. As another evidence, we have confirmed that 
triangular and sinusoidal shaking fields (that at the same 
fundamental frequency induce similar vortex velocities) 
produce similar results. 

The simulated system describes a very small region 
of the real macroscopic sample, where the mean vor- 
tex density and mean current density can be assumed 
to be homogeneous. In the experiment, the shaking ac 
field induces macroscopic currents J s h(t) that move vor- 
tices, that have been introduced in our simulations by 
a Lorentz force term F L = 4>qJ s Ji X z, where z is the 
versor parallel to the vortex direction. The explicit func- 
tional relation between J s h(t) and the shaking h ac (t) is 
not straightforward. In order to gain some intuition, we 
have applied the method developed in Ref. to obtain 
the distribution of J(t) for a certain geometry, assum- 
ing a non linear electric field-current constitutive relation 
E(J) = E C (J/ J c ) n —. We have chosen a simple geome- 
try (a disk) and we have estimated n by comparing non 
linear susceptibility measurements and polar plots of x" 
vs x' ■ Under this non linear regime, given an external 
h s h(t) = h s h cos(wi) the current throughout most of the 
sample can be described by a square waveform of fre- 
quency uj and normalized amplitude Jo(w)/J c > 1 (Fig. 
|3ta)), being the frequency dependence logarithmic. Thus 
as a first (and very crude) approximation we will simulate 
our model with J s h{t) = Jo(w)Sign(sin(wt)). 

In our simulations, we consider N v rigid vortices 
with coordinates in a two-dimensional rectangle of 
size L x x L y that evolve according to the dynam- 
ics rj-Vi = F l - V + F^ p + F L , where v, ; its velocity, 
r] the Bardeen-Stephen viscosity coefficient, FJ V = 

iij is the vortex- vortex interaction 
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FIG. 3. (a) Theoretical current density J/J c at a fixed po- 
sition, in a superconducting disk of radius a and a perpen- 
dicular sinusoidal magnetic field, as a function of t/P for two 
different frequencies, lcD, solid line and IOOcj, dashed line, 
w = f ■ Throughout most of the sample, the current can 
be described by a square waveform, (b) Calculated vortex 
mean velocity v, obtained from molecular dynamic simula- 
tions, as a function of t/P. 



and Ff = - J2k=i Fle- (Zl ^ ] fa-H*) is the pinning 
attraction. Here, </>o is the quantum of magnetic flux, A is 
the London penetration length, K i is the special Bcsscl 
function and f vv is a dimensionless parameter that can be 
related to the stiffness of the vortex lattice. The N p pin- 
ning centers are located at random positions R^, with 
strength and range F? and r p .We measure lengths in 

.2 

units of A, forces (per unit length) in units of /o = 87 ^° A3 , 
time in units of to = V^/ fo- We will consider N v — 1600, 
L x = 40A, L y = V3L x /2, N p = 25N V , r p = 0.2A and 
FE chosen from a Gaussian distribution with mean value 
F p = 0.2 and a standard deviation of 0.01F P . The equa- 
tions of motion are integrated using a standard Euler al- 
gorithm with step h = 0.04io, and a hard cut-off A = 4A 
for the vortex-vortex force. The mean vortex velocity 
is defined as v(t) = jpYli v V «W> where Vi(t) is the in- 
stantaneous velocity of the i- vortex in the force direction. 
The observable used to characterize the degree of order is 
the proportion of vortices with 6 neighbors P§ = 1 — rid, 
where is the density of disclinations. P§ is calculated 
at the end of each cycle. 

We assume that the main effect of temperature near 
the PE is the increase of the ratio f vp /fw, causing a 
spontaneous disordering of the VL. In the simulations, 
starting with a perfect ordered lattice (i.e. Pq = 1) and 
leading the system to evolve without any applied external 
force, there is a spontaneous creation of disclinations, and 
a decrease in P§ . In the inset of Fig. [U the spontaneous 
Pq starting from an ordered configuration is plotted as 
a function of f vp / f vv . We identify the parameters corre- 
sponding to the experimental temperature region shown 
in Fig. [1] as the range where the VL spontaneously dis- 
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FIG. 4. Pe versus t/P starting from an ordered and a disor- 
dered configuration, (a) For two different to and fixed /t. (b) 
For two different /l and fixed cj. (Inset) P§ versus fvp/fw, 
showing a spontaneous creation of disclinations, and a de- 
crease in Pe at around f vp /fw = 0.25. 



orders. Therefore, we set the ratio f vp /fw = 0.25 to 
simulate the experimental system at a fixed temperature 
(experiments of Fig. [2}. The dc depinning force per unit 
length F^ c — (f>oJ c ~ 0.01/ was obtained by increasing 
a dc Fl until the mean velocity v exceeds a "criterion 
voltage". The shaking frequencies have been selected in 
the low frequency regime (uj « u c ~ — ~ 0.11}^, 
where pinning forces dominate over viscous drags, be- 
cause the experimental shaking frequencies are well be- 
low f c = u c /2tt > IMHz 

Figure [3jb) shows the calculated mean velocity as a 
function of time during two shaking cycles for two differ- 
ent frequencies. After switching the force direction, v de- 
cays towards a steady regime characteristic of a dc driv- 
ing force, similar to that reported in Ref. Q that could 
be associated to a FSS. The transient time r and the final 
mean velocity depend on the shaking force. At very low 
frequencies, the system reaches the steady regime inside 
each cycle (black curves in Fig. [3jb)) , and the response 
is essentially dc. This is not the case at higher frequencies 
(gray curve), were we observe that both r and the final 
mean velocity slightly vary while the system evolves from 
a cycle to the next one. This fact is more clear in Fig. [4j 
that shows Pq as a function of the number of cycles of the 
shaking force, starting from different initial conditions, at 
a fixed amplitude and two different frequencies Qa)) and 
a fixed frequency and two amplitudes slightly above the 
depinning force ((4(b)). After a new transient time r oc , 
that can involve many shaking cycles, the disclination 
density remains fluctuating around a steady value, in a 
dynamic state independent of the initial condition, that 
depends on the amplitude of the shaking force. When 



the force is very near the dc depinning force, a clear fre- 
quency dependence is also observed (Fig. Ufa)). It be- 
comes clear that the higher the shaking frequency and 
amplitude of the shaking force (both effects expected by 
increasing, in the experiment, the frequency of a shaking 
field), the more ordered the resulting VLC, in agreement 
with that observed in our experiments. In this process, 
there are two different transient times that involve the 
creation or annihilation of disclinations: the characteris- 
tic dc time to decay to the FSS inside a cycle r ( prob- 
ably related with the experimental transient time t-i in 
Ref. [2j ) and the transient time r ac that can take many 
shaking cycles and could be related with the transient 
response observed in ac transport experiments 2,6 . This 
second process is inherent of an oscillatory drive, and dy- 
namically reorganizes the system in a VLC independent 
of the initial conditions. 

In conclusion, a dynamic reorganization in the plas- 
tic regime that keeps memory of the frequency drive has 
been observed experimentally and numerically. Through- 
out the order disorder crossing phase boundary in NbSe2 
there is a large collection of VLCs, that can be accessed 
after the application of ac shaking fields of different fre- 
quency, either ordering or disordering the initial state. 
The order of the VLCs is determined by non invasive ac 
susceptibility measurements, showing that higher shak- 
ing frequencies lead to more ordered VLCs. The system 
forgets the initial condition and but keeps memory of the 
shaking frequency, suggesting that the nature of the fi- 
nal states is more complex than previously conjectured--. 
Molecular dynamics simulations performed for a compat- 
ible set of interaction parameters in the plastic regime 
have reproduced qualitatively the salient features of ex- 
perimental shaking protocols, revealing a plausible na- 
ture of the experimental VLCs. The dynamic steady 
states are reminiscent of FSS observed in other systems, 
as VL disclinations and mean vortex velocity fluctuate 
around a steady value. As a novel fact, we have identi- 
fied two characteristic time constants: one characteristic 
of the dynamics inside each cycle (dc) and the second 
related to ac dynamics. Furthermore, as the theoreti- 
cal model is similar to the one used in the simulations 
of other interacting systems^, we expect that frequency 
dependent final states should be observed experimentally 
in a wide collection of particle systems with plastic flow 
regimes. A deeper study of the transient processes and its 
connection with the depinning transition will be object 
of future work. 
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